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Constraining the topology of the Universe using the 
polarised CMB maps 
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ABSTRACT 

We study the possibility for constraining the topology of the Universe by means of the 
matched circles statistic applied to polarised cosmic microwave background (CMB) 
anisotropy maps. The advantages of using the CMB polarisation maps in studies of 
the topology over simply analysing the temperature data as has been done to-date are 
clearly demonstrated. We test our algorithm to search for pairs of matched circles on 
simulated CMB maps for a universe with the topology of 3-torus. It is found that the 
noise levels of both Planck and next generation CMB experiments data are no longer 
prohibitive and should be low enough to enable the use of the polarisation maps for 
such studies. For such experiments the minimum radius of the back-to-back matched 
circles which can be detected are determined. We also showed that the polarisation 
generated after reionisation does not have an impact on detectability of the matched 
circles. 
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1 INTRODUCTION 

According to General Relativity the local properties of 
spacetime geometry are described by the Einstein gravi- 
tational field equations. However, they do not specify the 
global spatial geometry of the Universe, i.e. its topol- 
ogy. This can only be constrained by observations. The 
question of a possible multiply connecte d topology of 
our c osmic space was already raised by ISchwartzschildl 
(1900) b efore the formu lation of Genera l Relat ivity as 
well by Ide Sitter! \l9lH ) and iFriedmannl \l924 ) imme- 
diately following the formulation of the theory. How- 
ever, because of a lack of observational data probing the 
largest scales of the Universe, the standard cosmological 
model was tacitly assumed to comprise a universe pos- 
sessing a simply-connected topology. One of the first at- 
tempts to constraining the Universe's topolog y was under- 
taken utilising the distribution of quasars by iFang fc Satd 
(1985). However, the first measurements of the CMB 
anisotropy led to a sig n ifican t improvement in constraints 
( Stevens. Scott. & Silk! 19931: de Oliveira-Costa fc Smood 



zon of the observable Universe and therefore provide the 
best chance for detection of the signatures of multi- 
connectedness. 

Moreover, in the last decade the study of topology has 
also drawn considerable attention due to various anoma- 
lies observed on the largest angular scales in the WMAP 
data. Some of these ano malies, such as the suppression of 
the quadrupole moment (Be nnett et al.ll2003h and an align- 
ment betwee n the preferred axes of the quadrupole and 
the octopole (Ide Oliyeira-Costa et alll2004l ; ICopi et alj2004 



iBielewicz et alj|2005h . may be explained as effects caused by 
a multi-connected topology. 

Current constraints on topology were established 
by searching for the signatures of multi-connectedness 
CMB sky maps either in the s pherical harmonic 



doma in (|de Oliveira-Costa et alj 12004 : Kunz et al 
2008f) or directly in pixel s pace ijCornish et al. 
Aurich Lustig. fc Steinerll2005l, 120061 ; lThedl2006l : iKev et all 



1 19951 ; Ide Oliveira-Costa. Smoot. fc Starobinskv 19961) Such 
anisotropy maps survey scales comparable to the hori- 
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20071 ; IBielewicz fc Bandavl l20lJ v In the latter case, 
one of the most sensitive methods is the search for 
matched circles in the temperature anisotropy patterns 
|Cornish. Spergel. fc Starkmanl Il998h . By means of this 
method as applied to the WMAP maps, a large class of 
topologies has been ruled out. However, the method is 
not inherently limited to temperature anisotropy studies. 
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It can also be applied to the CMB polarisation data. In 
this work, we investigate such an application of the method. 
We test it on simulated CMB maps for a flat universe 
with the topology of a 3-torus, and explicitly consider the 
possibility for the detection of matching circle pairs for 
data with an angul ar resolution and noise leve l character- 
istic of the Planck JPlanck Collaboration ||2005h and COrE 
l|COrE Collaboration 1120111 ) data. The latter will be treated 
as a reference mission for the next generation of CMB ex- 
periments. 

The paper is organised as follows. In Section[2]we briefly 
review some basic properties of the CMB temperature and 
polarisation anisotropies. A description of the simulations 
for a universe with multi-connected topology used to test 
the reliability of our codes is presented in Section [3] Sec- 
tion U is devoted to a description of the statistic adopted in 
our studies. The disadvantages of the search of the matched 
circles in temperature maps are presented in Section [5] and 
in Section [6] the application of the method to simulated po- 
larisation maps is shown. In the last section we present our 
conclusions. 

2 CMB ANISOTROPIES 

In this section we provide a brief resume on the properties 
of CMB temperature and polarisation anisotropy. 

2.1 Temperature anisotropies 

The CMB is observed as blac kbody radiation wit h tem- 
perature T = 2.725 ± 0.002 K jMather et al.lll999T ). After 
accounting for the peculiar motion of the Earth with re- 
spect to the CMB rest frame, the radiation is seen to ex- 
hibit tiny temperature fluctuations of order AT /To ~ 10 -5 . 
In the Newtonian gauge, the temperature fluctuations in a 
given direction of the sky, n. can be written as 

^ (n) = ($L + n-v lTr + £ e~^\ip + 4>)dr , 

(1) 

where <5 7 is the fractional energy density fluctuations in the 
CMB, v is the baryon velocity, cf> and tp are the two gravi- 
tational potentials in terms of which the metric fluctuations 
are given by ds 2 = a 2 (r) [-(1 + 2ip)dr 2 + (1 - 2<f>)dx 2 } , and 
the optical depth «(r) = a(r')aTn e x e dr' . r is conformal 
time with subscripts and r denoting the present time and 
the time of recombination, respectively. 

The different terms in equation (JTJ) describe different 
physical effects. The first term corresponds to photon en- 
ergy density fluctuations and gravitational potential redshift 
(so called the Sachs- Wolfe effect, SW), the second to the 
Doppler effect caused by baryon motions at recombination 
and the third is the so-called the integrated Sachs- Wolfe 
effect (ISW) caused by the variation of the gravitational po- 
tential with time. The equation does not include the terms 
coming from other secondary effects such as the Sunayev- 
Zeldovich and gravitational lensing effects. However, they 
are important for angular scales smaller than the scales of 
interests in this paper (~ 20'). 



2.2 Polarisation anisotropies 

Linear polarisation in the CMB is generated by Thomson 
scattering of photons by electrons either at the moment of 
recombination, when photons decouple from the primordial 
plasma, or during reionisation when, due to partial ioni- 
sation of the matter, free electrons reappear again in the 
Universe. The former produces a pattern of polarisation 
anisotropy at small angular scales while the latter generates 
structure on large scales, corresponding to the angular scale 
of the Universe at the moment of scattering. The small-scale 
polarisation anisotropy is considered to be a snapshot of the 
last scattering surface (LSS). 

To generate polarisation the radiation incident on the 
electrons needs to have a quadrupole moment. It can be 
produced either by scalar or tensor perturbations. Depend- 
ing on this, the polarisation pattern of the sky will mani- 
fest different types of parity. The polarisation modes gen- 
erated by the former will show "electric-type" parity while 
the modes generated by the latter - "magnetic-type" parity 
<|Zaldarriaga fc Seliak|[l997l ). As is common, we will refer to 
these two types of polarisation components as the E-modes 
and B-modes, respectively. We will not consider in this paper 
the impact of the gravitational lensing effect on the primor- 
dial polarisation since it is not important for the angular 
scales of interest here. 

In the case of scalar perturbations, gradients in the ve- 
locity field, diVj, of the photon-baryon fluid on the LSS 
lead to a quadrupole component of the temperature seen by 
the scatterer AT(n') « n'ifi'jdiVj , which, through Thomson 
scattering, is converted into polarisation. Then, the polari- 
sation expr essed in terms of the Stokes parameters Q and U 
is given by l|Seliak fc Zaldarriagalll99§ ) 

(Q+iU)(n) oc ot j dn (m-n) 2 AT(n) oc AT r rn m? diV j\ Tr , 

(2) 

where ot is the Thomson cross section, Ar r is the width of 
the LSS, m = ei + ie2 and ei, e.2 are directions perpen- 
dicular to h that are used to define the Stokes parameters. 
Tensor m % m? assures that the quantities Q ± iU transform 
under rotations of (ei, e.2) by an angle 7 as spin 2 quantities 
i.e. (Q ± it/)'(n) = e* 2i ~t(Q ± iU)(n) . 

The Stokes parameters are direct measurables. How- 
ever, because of their spin 2 nature, they are not conve- 
nient for the analysis of the polarisation field, for which the 
use of spin maps of the E-modes and B-modes is pre- 
ferred and allows easy discrimination of the type of pertur- 
bations generating the polarisation anisotropy. Such quan- 
tities can be related to the Stokes parameters by means 
of spin-raising 9 and spin-l owering 6 differential operators 
|Zaldarriaga fc Seliaklll997r i. i.e. 

E(ft) = ~[m(Q+iU)(fi)+m(Q-iU)(n)] , 

B(n) = L[BB(Q + iU)(n)-3d(Q-iU)(n)] . (3) 

It is worth noticing that the former equation corresponds to 
the divergence operator, the latter to the rotation operator 
of the spin-2 fields on the sphere. Thus, these relations can 
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be seen as analogs of Gauss's law and Ampere's lawQ re- 
spectively, in electrodynamics, and the E-mode and B-mode 
maps as analogs of the electric charges and the current den- 
sities. This explains why these two components of the polar- 
isation field are often called the "electric" and "magnetic" 
modes. 

To avoid the direct differentiation of the Q ± iU maps 
decomposition can proceed in the spin 2 spherical harmon- 
ics basis ±iYt m i.e. (Q±iU)(n) = Y^t, m a ±2,«m ±2^m(n). 
Then, the relations ([3J expressed by the appropriate spher- 
ical harmonics coefficients take the form 

E _ "2,(m + 1-2,(m B _ ■ a 2,im ~ 1- 2,£m i A \ 

a lm — g ' aira ~ 1 2 ' 

The af m and af m coefficients are then related to the E-mode 
and B-mode maps by 

E{h) = YlJ fet. 2) T a * m Yirn ^ ' 

Bin) = Eyf^|| a ?- Y ^(") • ( 5 ) 

It is important to note that the af^ 3 coefficients are not 
exactly equal to the spherical harmonic decomposition co- 
efficients of the E and B-modes maps. The use of the af^ 3 
coefficients without the factor \fil + 2)l/(£ — 2)! is partly 
a matter of convention. It is motivated by the fact that 
the power spectrum computed from such coefficients is the 
same on small scales as for the Stokes parameters them- 
selves, while that of the E and B maps differs by a factor 
of oc £ 4 . One consequence of the yjjh + 2)1/ (I — 2)! factor is 
that white noise for the Stokes parameters anisotropy fields 
will transform to coloured noise for the E and B maps. 



In the simulations, we did not take into account the 
B-modes for multi-connected universes. The amplitude of 
the polarisation generated by gravitational waves is much 
smaller than the E-mode polarisation and the noise level for 
the Planck data is too high to allow the use of these maps 
for constraining the topology. We also do not consider effects 
that generate B-mode polarisation from the gravitational 
lensing of the E-modes. As is the case for any secondary ef- 
fect, this dilutes signatures coming from the topology. Nev- 
ertheless, it is important for scales much smaller than those 
considered in this work (~ 20'). However, it is important to 
note that in the simulated maps used in this paper B-modes 
are present due to the noise. 

To study the signatures of a given topology, a CMB 
map is required with resolution comparable to the angular 
size of the first significant acoustic peaks of the polarisation 
map i.e. ^ max ~ 400 — 500. We have adopted £ max = 500 
in our simulations. The dimension of the fundamental do- 
main of the 3-torus was L — 2 c/ Ho, which is about three 
times less than both the diameter of the observable Universe 
i.e. ~ 6.6 c/ Ho, and current lower bounds on the size of the 
Universe obtained fro m a matched circles study of tempera- 
ture anisotropy maps (|Kev et alj|2007l ; iBielewicz fc Bandavl 
12011ft . Although such a small 3-torus does not fit very well 
to the data, the simulations are used only to demonstrate 
and compare the performance of the statistic for matched 
circles over a wide range of radii. A universe where there 
are many pairs of matched circles of different radii serves 
very well for these objectives. The values of the cosmologi- 
cal parameters adopted corresponded to the A.CDM model 
(see lLarson et al1l201ll . Table 3) determined from the 7-year 
WMAP results. The time needed for the generation of one 
such CMB map on a single processor with clock speed 3 GHz 
is about 42 hours. 



3 SIMULATIONS OF THE CMB POLARISED 
ANISOTROPY MAPS IN THE 
MULTI-CONNECTED UNIVERSE 

To test the reliability of the codes used to search for the 
signature of a multi-connected topology in the Universe, 
we performed simulations of CM B skies for a flat un iverse 
with the topology of a 3-torus (|Riazuelo et al.l 12004). The 
first step in this process involved the computation of the 
coefficients of the E-modes, af m , in an analogous manner 
to the temperature anisot ropy coefficients as described in 
IBielewicz fc Bandavl (|201lft . The only significant difference 
was the use of the response function for the E-mode po- 
larisation, Af (k, At), instead of that for the temperature 
in their equation (1). The coefficients were then convolved 
with an appropriate beam profile, and subsequently trans- 
formed to the experimentally observed Stokes parameters 
using the relation @ and inverse of the spin 2 spherical 
harmonics transform described in the previous section. Fi- 
nally, homogeneous noise was added to the Q and U maps 
at a characteristic level for a given experiment. 



1 because the measured Stokes parameters are independent of 
time we do not have time dependent terms in this case 



4 STATISTIC FOR THE MATCHED CIRCLES 

If light had sufficient time to cross the fundamental cell, an 
observer would see multiple copies of a single astronomical 
object. To have the best chance of seeing 'around the uni- 
verse' we should look for multiple images of distant objects. 
Searching for multiple images of the last scattering surface 
is then a powerful way to constrain topology. Because the 
surface of last scattering is a sphere centred on the observer, 
each copy of the observer will come with a copy of the last 
scattering surface, and if the copies are separated by a dis- 
tance less than the diameter of the last scattering surface, 
then they will intersect along circles. These are visible by 
both copies of the observer, but from opposite sides. The two 
copies are really one observer so if space is sufficiently small, 
the CMB radiation from the last scattering surface will con- 
tain a pattern of hot and cold spots that match around the 
circles. 

The idea of using such circles to stu dy topology is due 
to ICornish. Spergel. fc StarkmarJ (|l99of ). Therein, a statis- 
tical tool was developed to detect correlated circles in all 
sky maps of the CMB anisotropy - the circle comparison 
statistic 



S Pt ria, <f),) 



(2X p (±0)X r (0 + 0„)) 
(X v {^) 2 + XriW) 



(6) 
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where {) = d<f> and X p {±4>), X r ((f> + <j>*) are temperature 
(or polarisation) fluctuations around two circles of angular 
radius a centered at different points, p and r, on the sky 
with relative phase </>*. The sign ± depends on whether the 
points along both circles are ordered in a clockwise direction 
(phased, sign +) or alternately whether along one of the cir- 
cles the points are ordered in an anti-clockwise direction 
(anti-phased, sign — ). This allows the detection of both ori- 
entable and non-orientable topologies. For orientable topolo- 
gies the matched circles have anti-phased correlations while 
for non-orientable topologies they have a mixture of anti- 
phased and phased correlations. The statistic has a range 
over the interval [—1, 1]. Circles that are perfectly corre- 
lated or anticorrelated have S = 1 or S = — 1, respectively, 
while uncorrelated circles will have a mean value of S = 0. 
To find correlated circles for each radius a, the maximum 
value 5'm a x( a ) = max S , ^ r .(a,0„) is determined. In case of 

anticorrelated circles the maximum value of —S^ r (a, 0») is 
used. 

In the original paper bv lCornish. Spergel. fc Starkmanl 

(1998) the above statistic was applied exclusively to temper- 
ature anisotropy maps. However, it can also be applied to 
polarisation data and in this work, we focus on its applica- 
tion to the E-mode map. In this case, the X map is simply 
the map of the E-mode. 

In order to speed up the computations, one can 
use the fast Fourier transform (FFT) along the cir- 



cles, x p (4>) 



Y, m X p , m exp(im<f>), and, by rewriting 



the statistic as Sp tT .(a, <j),) = J? s m exp(— im0»), where 
s™, = 2E ra ^,>7.^,»/E„ {\ x p^\ 2 + \X r , n \ 2 ), use the in- 
verse fast Fourier transform. It reduces the computational 
cost 5*jJ; r (a, (/>,) from ~ N to ~ 7V 1//2 log N operations, where 
N is the number of pixels in the map. 

In our studies w e will use the statistic as modified in 
ICornish et all l|2004l ) to include weighting of the mth har- 
monic of the temperature around the pth circle X v>ra by the 
factor |m| such that the number of degrees of freedom per 
mode are taken into account. 



Sp, r (a, 4>*) = 



2E 



I | , m ^r, m & 



En M (l*p,. 
where the Fourier coefficients X„ 
in the case of the S~ 



+ \Xr, n \ 2 ) 



(7) 



t P)m are complex conjugated 
■ „ r statistic. Such weighting enhances 
the contribution of small scale structure relative to the 
large scales fluctuations. It is especially important because 
the large scale fluctuations, which are dominated by the 
integrated Sachs-Wolfe (ISW) effect for the temperature 
anisotropy and the reionisation effect in the case of polarisa- 
tion, can obscure the image of the last scattering surface and 
reduce the ability to recognise possible matched patterns on 
it. 

The general search explores a six dimensional space: 
the location of the first circle centre, (9 P , (j> p ), the location of 
the second circle centre, (8 r ,(f> r ), the angular radius of the 
circle a, and the relative phase of the two circles </>„. The 
number of operations needed to search for circle pairs scales 
as ~ TV 3 log A''; thus implementation of a matched circles 
test is computationally very intensive. As a consequence, 
we restrict our analysis to a search for pairs of circles cen- 
tered around antipodal points, so-called back-to-back circles. 
Then, the search space can be reduced to four dimensions 
and the number of operations to ~ N 2 log TV. Although the 



back-to-back search can detect a large class of the topologies 
that we might hope to find, there remain those topologies 
that predict matched circles without antipodal configura- 
tions such as the Hantzsche-Wendt space for flat universes 
or th e Picard space for hyperbolic universes |Aurich et alj 

12004) . 

As in lBielewicz fc Ban day (201l]), we used the healpijsQ 
scheme with a resolution parameter N s id e = 512 to define 
our search grid on the sky. To accommodate the use of the 
FFT approach, n = 2 r+1 points are used for each circle, 
where r corresponds to the resolution parameter of the map 
given by N si d e = 2 r . By definition this is also the angular 
resolution used to step through </),. For a, we used steps a 
little bit smaller than the characteristic scale C of coherent 
structures in the map, i.e. 2 c /3. As we discuss below, to 
increase the signal-to-noise ratio the polarisation maps were 
smoothed with a gaussian beam profile with FWHM of 50'. 
Then, the scale of coherence is defined relative to the two- 
point correlation function as twice its angular width at half 
the maximum of the zero lag value. For polarisation maps 
this is around 20'. Contrary to the temperature anisotropy 
this scale is smaller than the smoothing scale due to the 
bluer power spectrum, dominated by small angular scales, 
of the polarisation maps. The values of the E-mode map 
at each point along the circle were interpolated based on 
values for the four nearest pixels with weights inversely pro- 
portional to the distances between a given point and the 
pixel centers. Other methods of interpolation can be used 
but we have verified that even using the value of the nearest 
pixel does not change the value of the statistic significantly. 
The time needed for a search of one map using this scheme 
on a single processor with clock speed 3 GHz is about 180 
hours. 

To draw any conclusions from an analysis based on the 
statistic 5'max( Q ) if is important to correctly estimate the 
threshold for a statistically significant match of the circle 
pairs. The chance of random matches is inversely propor- 
tional to the number of coherent structures along the cir- 
cles, therefore a false positive signal level of the statistic is 
especially large for circles with smaller radii. We used sim- 
ulations of the maps with the sam e noise properties and 
smoo thing scales as the Planck data (|Planck Collaboration I 
2005) to establish the threshold such that fewer than 1 in 
100 simulations would yield a false event. 

It is important to note that the false detection level 
is smaller for higher resolution maps. Co nversely, as shown 
by ICornish. Spergel. fc Starkmanl (|l998T ). the value of the 
statistic for matched circles 



(8) 



depends on the ratio of the signal rms cr s to the noise rms cr n 
ratio, denoted £ for the map, thus the efficiency of the statis- 
tic to detect matched circles is increased by smoothing the 
data. A smoothing scale should therefore be adopted that is 
a reasonable trade-off between these two requirements. 

To eliminate the regions most contaminated by Galac- 
tic foreground residuals, we utilise the p olarisation analysi s 
mask P06 defined by the WMAP team l|Page et all 120071 ). 



http: / /healpix. jpl.nasa.gov 
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However, the statistic is very sensitive to the masking, par- 
ticularly if a significant fraction of one or both circle pairs 
lies in the masked region. In this case, there is a significant 
probability to find a chance correlation of the temperature 
pattern between unmasked parts of the circles, thus increas- 
ing the false detection level. The effect is most pronounced 
for circles with the largest radii, close to 90 degrees, as well 
as for very small radii. To avoid this, we restrict our analysis 
to those pairs for which less than half the length of each cir- 
cle is masked. Although the statistic computed in this way is 
not optimal, it is much more robust with respect to masking. 



5 DISCUSSION OF DEGRADING EFFECTS 
FOR TEMPERATURE MAPS 

Hitherto, the search for matched circles was performed only 
for the CMB temperature anisotropy maps. The main ob- 
stacle to the the application of this method to polarisation 
data is the very high noise level of the currently available 
CMB polarisation maps. This situation will change with the 
release of the Planck data. As we will show below, in this 
case the typical noise amplitude will be sufficiently low to 
enable the use of the polarisation data in studies of topology. 
This provides an opportunity to avoid some of the problems 
inherent in the analysis of temperature maps. 

Since the signatures of topology are imprinted on the 
surface of last scattering, any effects that dilute this im- 
age will also degrade the ability to detect such signatures 
by means of the matched circles statistic. In the case of the 
temperature fluctuations, there are two sources of anisotropy 
generated at the LSS: the combination of the internal pho- 
ton density fluctuations and the SW effect, described by the 
first term of equation Q, and the Doppler effect, described 
by the second term. In the latter case, the correlations for 
pairs of matched circles can be negative in universes with 
multi-connected topology. As we can see in Fig. [TJ for pairs 
of back-to-back circles with a radius smaller than 45° the dot 
products n v have opposite signs for two nearly antipodal di- 
rections of observations. The Doppler term becomes increas- 
ingly anticorrelated for circles with smaller radius. One can 
show that on average in a flat universe with orientable topol- 
ogy the statisti c © for this term is (S~(a)) = — cos(2a) 
(Riazuclo ct al. 2006). Anticorrelation of the Doppler term 
significantly decreases correlations with the circle radius for 
the total CMB anisotropy map. We can see this effect in 
Fig-0f° r a simulated CMB sky in a universe with the topol- 
ogy of the 3-torus. 

The ISW effect, as for any other secondary anisotropy 
contribution, will also smear out signal coming from the LSS. 
However, because it dominates only at large angular scales, 
this can be remedied either by high-pass filtering of the map 
or by using the matched circles statistic with m weighting 
to prevent the statistic from being dominated by the 
large angular scales. 

The consequences of these degrading effects are weaker 
constraints on the topology of the Universe obtained from 
the matched circle statistic. As we can see in Fig. [2] use of 
the CMB map without both of the degrading effects would 
allow us to impose lower bounds on the minimum radius of 
the correlated circles which can be detected much below th e 
present constraints a m i n ~ 10° IjBielewicz fc BandavH201lT ). 



SEARCH OF MATCHED CIRCLES IN 
POLARISATION MAPS 



As already mentioned in Section \2. 2 1 a polarisation map at 
small angular scales can be considered as a snapshot of the 
LSS. Theoretically, then, the polarisation provides a better 
opportunity for the detection of multi-connected topology 
signatures than a temperature anisotropy map. The only 
serious issues preventing its use in studies of topology are in- 
strumental noise and the correction of the polarised data for 
the Galactic foreground. Since the polarisation anisotropy 
amplitude is around ten times smaller than the amplitude for 
temperature, these constitute more serious problems than 
for temperature anisotropy maps. 

In addition, since polarisation is generated by the scat- 
tering of a quadrupolar distribution of radiation on free elec- 
trons at the moment of recombination or during reionisation, 
the orientation of the local quadrupole can be a significant 
problem to this analysis. The orientation of the dipole and 
quadrupole are determined uniquely by the initial density 
field. Thus, one can see from equation © that the generated 
polarisation depends on the projection of the quadrupole 
on to the plane perpendic ular to the line of sight. As shown 
by iRiazuelo et al.l {2006) this projection can be expressed 
by 1 — (n • v) 2 . This is analogous to the projection of the 
dipole direction on the line of sight for the Doppler term, as 
discussed in the previous section. Clearly, the contribution 
is included in our simulations, and any deleterious effect is 
accounted for in the statistical analysis. 

The other effects that might degrade the signal from the 
LSS are the polarisation generated after reionisation and the 
gravitational lensing effect. The former prevails at very large 
angular scales where, as for the ISW effect and temperature 
anisotropy data, its influence can be eliminated by either 
high-pass filtering or the use of the statistic 0. We will 
investigate this below in subsection 16. II The lensing effect is 
most significant at very small scales (~ 3') and is negligible 
for the range of scales of interests in this paper. 

We have tested the perfor mance of the matched circles 
statistic for simulated Planck dPlanck Collaboration "1|2005T ) 
and COrE {COrE Collaboration I 2011 ) polarisation maps. 
Although the latter has not currently been selected as a 
candidate mission by ESA, we will consider it as a useful 
reference mission for the next generation of observational 
CMB projects. 

Before we present the results for simulations, it is worth 
considering a rough estimation of the performance of the 
matched circle statistic for different experiments. As we have 
already mentioned, the noise level of the WMAP polarisa- 
tion data is too high to enable detection of the matched 
circles. Even for the map assembled from the WMAP Q, V 
and W-band polarised data with inverse noise weighting we 
find Smax ~ 0.05. This is considerably below the false de- 
tection threshold of the statistic. There is a much improved 
prospect for the detection of the signatures of topology for 
the Planck data. As for the WMAP data, to decrease the 
noise level we consider the map obtained from a combination 
of the data with the lowest noise levels using inverse noise 
variance weighting. Using the sensitivity and resolution of 
the polarised anisotropy map derived from the 100, 143 and 
217 GHz frequency maps, i.e. with A%' u 59 /iKar cmin 
and #fwhm = 10 arcmin l|Planck Collaboration II2005T ). re- 
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Figure 1. The baryon velocity v at the last scattering surface seen from two different directions n, n' in the multi-connected universe. 
For the radius of the matched circles larger than 45° the projections of the baryon velocity on the direction of observations (dashed 
lines) for the observers O and Oi are strongly correlated. For the radius equal to 45° (observers O and O2) the velocity projections are 
independent and for the radius smaller than 45° (observers O and O3) the projections are anti-correlated. The projections are denoted 
by blue and green arrows. 




Figure 2. An example of the S^ax statistic for a simulated CMB temperature anisotropy map of universe with the topology of a cubic 
3-torus with dimensions L = 2 c/Hq. The dotted, solid, dashed and three dot-dashed lines show the statistic for CMB maps of the 
Sachs- Wolfe (SW) effect, the positive and negative correlations of the Doppler effect and total anisotropy, respectively. 



spectively, we get S max ~ 0.6. This is sufficient to enable the 
detection of pairs of matched circles. Conversely, the noise 
level for single frequency maps seems to be too large to allow 
such a test - even the 1 43 GHz channel with a sens itivity 
A%' u « 82 ^Karcmin (jPlanck Collaboration 1 I2005T I only 
renders a value of Smax ~ 0.4. The best chance of detec- 
tion for such signatures of topology arise for an experiment 



with specifications similar to those of the COrE project. 
For a channel with minimal contamination by the Galac- 
tic foreground, i.e. 105 GH z, with a polarisation sensi tivity 
of A%' u » 4.6 ^Karcmin (|COrE Collaboration I lioTlh , the 
statistic Smax ~ 0.99 achieves a value close to the maximum. 

As described in Section[3] we simulated the CMB polar- 
isation maps, in terms of the observable Stokes parameters 
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Q and U, for a universe with the topology of a 3-torus of 
size L — 2 c/ Hq. The resolution of the maps and the noise 
rms correspond to the Planck 100, 143 and 217 GHz fre- 
quency maps coadded with inverse noise weighting and the 
COrE 105 GHz frequency map. In both cases, the maps were 
convolved with a gaussian beam of FWHM=10' and homo- 
geneous noise with a sensitivity of A%' u ~ 59 /^Karcmin 
and A^' 17 w 4.6 /iKarcmin, respectively, were added. 

Since the direct use of the Stokes parameters in the 
analysis would require the rotation of values into a com- 
mon coordinate system (defined by the pairs of points in 
question) which can significantly slow down the computa- 
tions, the analysis was carried out for E-mode maps ob- 
tained from the Q, U maps by means of the relation Q. 
The use of the spin maps allows us to speed up compu- 
tations by the application of an FFT along the circles, as 
described in Section [4] The maps were also smoothed with 
a beam of FWHM=50' to increase the signal to noise ra- 
tio. Because of the very blue power spectrum of the E-mode 
maps, caused by the + 2)!/(^ — 2)! factor (see equation 
©), in order to remove the noise on the smallest angular 
scales the width of the beam profile has to be wider than 
would be the case for the temperature map. The smoothing 
serves also to eliminate from the analysis multipoles of order 
larger than the maximum order used for the simulated CMB 
maps i.e. ^ max = 500. 

The S^ax statistic for the simulated maps is shown in 
Fig. [3] As expected the amplitudes of the peaks do not de- 
crease with the radius of the circles as in the case of the 
temperature anisotropy maps. We see that pairs of matched 
circles can be detected for the Planck coadded map. How- 
ever, comparing with Fig. [2] one should notice that the rel- 
ative amplitude of the peaks with respect to the average 
correlation level for the circles with small radius is not big- 
ger than for the temperature map. Thus, the constraints on 
topology will not be much tighter than those derived from 
an analysis of temperature maps. A much better perspective 
arises for the COrE maps. The signal of the multi-connected 
topologies is very pronounced in this case enabling the de- 
tection of matched circles with very small radius. It is also 
worth noticing that the width of the peaks is smaller than 
for temperature maps. It is a consequence of domination of 
the small angular scales in the power spectrum of the po- 
larised maps leading to more narrow peak of zero lag of the 
two-point correlation function. 

To determine more precisely the minimum radius of 
the detectable circles it is needed to estimate the false de- 
tection level of the statistic. We used for this purpose 100 
simulations of the Planck coadded map for a universe with 
a simply-connected topology. The detection level was esti- 
mated for both the full sky Q and U maps and the masked 
maps. The former corresponds to the case of a perfect sepa- 
ration of the CMB and the Galactic foreground polarisation 
and yields optimal performance of the statistic. The latter 
presents a more realistic estimation taking into account the 
necessity of removing those regions most contaminated by 
the Galactic foregrounds. 

We used m asks based on th e P06 mask derived by the 
WMAP team l|Page et al.1 l2007h for their analysis of po- 
larised data. In order to minimize ambiguous modes that 
arise from the decomposition of the Stokes parameters into 
the E and B-mode on incomplete sky coverage, we used the 



following procedure for computation of the E-mode maps. 
The P06 mask was extended an d apodized employing the 
method described by iKiml (|201ll ) (referred to as a 'process- 
ing mask' therein) i.e. the mask was widened so that after 
apodization with a Gaussian profile it is zero in regions de- 
fined as such by the original mask up to a precision of 10 -6 . 
Then, the E/B decomposition of the Q and U maps was 
performed with this mask. As for a full sky analysis, in or- 
der to increase the signal to noise ratio, the E-mode maps 
were also smoothed with a beam of FWHM=50'. Due to 
the apodization, most of the ambiguous modes are localized 
close to the mask edges. They can then be eliminated from 
the computation of the statistic by use of an additional mask 
applied to the derived E-mode map, but further broadened 
then apodized in a similar fashion to the previous mask. In 
this case, apodization helps to minimize the ringing effect 
which might appear on account of the use of the FFT in 
computations of the matched circles statistic. In both cases 
we used apodization with a Gaussian of FWHM=60'. Then, 
the fraction of the unmasked part of the sky for the cut em- 
ployed in the computation of the statistic is f s ^ y = 0.55. 
It is much less than for the original mask with the fraction 
/sky = 0.73. The masks are shown in Fig. [4] Such an aggres- 
sive excision of data corresponds to a rather extreme case. 
We may hope that the broader frequency range and lower 
noise level of the Planck data will enable us to significantly 
extend the useful sky coverage for cosmological analysis. 

This procedure for the computation of the E-mode sig- 
nal does not remove completely leakage from the B-modes. 
However, the residual is no larger than the noise level in 
smoothed polarisation maps that is negligible compare to 
the prevailing of E-mode signal. It is worth noticing that in 
our procedure the residual signal from B-modes is addition- 
ally damped by smoothing. 

The false detection levels established from the require- 
ment that fewer than 1 in 100 simulations should yield a 
false match are shown in Fig. [3] corresponding to two ex- 
treme cases: analysis on the full sky and on a map with a 
very conservative mask. The detection threshold for both 
Planck and COrE should lie somewhere between these two 
curves. However, more detailed studies of the influence of the 
Galactic foreground polarisation is required, but is beyond 
the scope of this paper. 

The minimum radius, a m i n , of the correlated circles 
which can be detected for a given map is determined by 
the intersection of the peaks in the matching statistic with 
the false detection level. The height of the peaks indicates 
that the minimum radius is around 20° for the Planck and 
around 5° for the COrE data in the case of the masked map. 
If we do not detect any pairs of matched circles, this radius 
will determine a lower bound on the size of the fundamental 
domain. 

6.1 Impact of reionisation on the matched circles 
statistic 

One of the secondary anisotropy effects which might degrade 
our ability to detect pairs of matched circles in polarisation 
induced by the Thomson scattering of the CMB photons 
on the reionised matter. To evaluate the level of impact of 
the effect on the performance our statistic, we simulated 
the CMB map for the multi-connected universe both with 
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Figure 3. An example of the Smax statistic for a simulated polarised CMB map of the universe with the topology of a cubic 3-torus of 
dimension L = 2 c/Ho- The solid and dashed lines show the statistic for simulated polarisation maps with angular resolution and noise 
level corresponding to the Planck and COrE data, respectively. The dot-dashed and three dot-dashed lines show the false detection levels 
for the statistic estimated from fOO Monte Carlo simulations of the Planck coadded fOO, 143 and 217 GHz frequency polarisation maps 
for the full sky and cut sky analysis, respectively. 




Figure 4. Masks used in the analysis. The dark blue region corresponds to the part of the sky unused in polarisation analysis, defined 
by the WMAP team in the P06 mask. The light blue and yellow bands correspond to the apodized regions of the extended mask used 
for the E/B decomposition and the mask used in the search for the matched circles, respectively. The red region corresponds to the 
unmasked and unsmoothcd part of the sky. 



and without reionisation and compared the statistics. The 
optical depth for the reionisation effect was set as k — 0.088. 
This comparison showed that the relative difference between 
the statistics is of order of 10~ 3 . Therefore, the influence 
of polarisation generated after reionisation is negligible for 
the detectability of the matched circles. This should not be 
a surprising conclusion since the S max statistic is sensitive 
rather to small angular scales while the polarisation effect 
of reionisation prevails at large scales. In case of the E-mode 
map the small scales are also additionally amplified by the 



\J(l + 2)\/(£ — 2)! factor, thus the large angular scales are 
not so significant as for the temperature anisotropy. 



7 CONCLUSIONS 

We have studied the possibility of using CMB polarisation 
maps for constraining the topology of the Universe by means 
of the matched circles statistic. We have shown that the 
CMB polarisation maps are free from the degrading effect 
of the Doppler term that contributes to the temperature 
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anisotropy and that consequently allows matched circles to 
be detected with a radius much smaller than for the CMB 
temperature maps thus providing tighter constraints on the 
topology. Moreover, the detection is robust with respect to 
the secondary CMB polarisation generated after reionisa- 
tion. However, the application of polarisation data to stud- 
ies of the topology of the Universe depends to a large extent 
on systematic effects such as noise and Galactic foreground 
residuals. 

Using simulations of the Planck polarisation data we 
demonstrated that the level of the instrumental noise should 
be sufficiently low to enable the detection of the signatures 
of a multi-connected topology, though it will not be low 
enough to improve constraints on the topology derived from 
the temperature maps. Nevertheless, it can serve as an im- 
portant cross-check of the latter. The best chance for the 
detection of the matched circles is afforded by a map ob- 
tained from the inverse noise weighted combination of the 
three least noisy Planck frequency maps i.e. 100, 143 and 
217 GHz maps. The polarised Galactic emission should also 
be the weakest for the first two of these channels enabling 
the best estimation of the CMB polarisation. 

All the advantages of the polarisation analysis can be 
fully exploited in the case of the planned next generation of 
full sky low-noise experiments. To demonstrate this, we used 
simulations with the resolution and sensitivity correspond- 
ing to the COrE mission. In this case the smallest radius of 
the detectable back-to-back matched circles in the polarisa- 
tion maps is significantly smaller than for the temperature 
maps. An estimation based on our simulations shows that 
the radius is around 5°. 

The detectability of the matched circles depends heav- 
ily on our ability to reliably separate the CMB polarisation 
from the polarised emission of the Galaxy. We restricted 
our studies on this issue only to the application of a mask 
intended to remove those areas of the sky most contam- 
inated by the Galactic foreground emission. We assumed 
that the residuals of the Galactic emission outside the mask 
can be sufficiently well modelled and subtracted from the 
map. However, it is not certain to what extent this will be 
possible. As is well known, removal of the polarised Galactic 
emission is a more challenging problem than for the temper- 
ature maps since there is no frequency band where the CMB 
polarisation dominates over the polarisation of the Galactic 
foreground even for a small fraction of the sky. 

We also have to bear in mind that a search for matched 
circles on a masked map is able to constrain only those uni- 
verses with such dimensions and orientation of the funda- 
mental domain with respect to the mask that allow the de- 
tection of pairs of matched circles. Since the probability of 
overlooking circle pairs may be significant for a mask that 
removes a large fraction of the sky, employing a cut as small 
as possible is one of the crucial issues. The analysis of a po- 
larised version of an internal linear combination-type (ILC) 
map on the full sky is also unlikel y to be a credib l e solu tion 
of this problem. As shown in lBielewicz fc Bandavl (|201 if ) for 
the temperature case, Galactic foreground residuals in the 
Galactic plane of an ILC map can cause significant spurious 
correlations between the circles. In the case of polarisation 
maps it will probably be an even more pronounced effect. 
However, it is to be hoped that the high quality polarisation 
data from the Planck satellite will enable a better under- 



standing of the Galactic polarised emission and therefore 
the use of a larger fraction of the sky for cosmological anal- 
ysis than for the WMAP data. 
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